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We apply an ab-initio approach to the nuclear structure of odd-mass nuclei straddling 48 Ca. 
Starting with the NN interaction, that fits two-body scattering and bound state data we evaluate the 
nuclear properties of A — 47 and A = 49 nuclei in a no-core approach. Due to model space limitations 
^vq , and the absence of 3-body interactions, we incorporate phenomenological terms determined by fits 

to A = 48 nuclei in a previous effort. Our modified Hamiltonian produces reasonable spectra for 
these odd mass nuclei. In addition to the differences in single-particle basis states, the absence 
of a single-particle Hamiltonian in our no-core approach obscures direct comparisons with valence 
effective NN interactions. Nevertheless, we compare the fp-shell matrix elements of our initial and 
modified Hamiltonians in the harmonic oscillator basis with a recent model fp-shell interaction, the 
GXPF1 interaction of Honma, Otsuka, Brown and Mizusaki. Notable differences emerge from these 
comparisons. In particular, our diagonal two-body T = matrix elements are, on average, about 
■ 800-900keV more attractive. Furthermore, while our initial and modified NN Hamiltonian fp-shell 

' matrix elements are strongly correlated, there is much less correlation with the GXPF1 matrix 

£Sj , elements. 
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I. INTRODUCTION 



The low-lying levels of the A=47-49 nuclei have long been of experimental and theoretical interest. On the one hand, 
extensive experimental information about these nuclei is available [H-Q an( A, on the other hand, this is a suitable 
nuclear mass region for developing and testing effective fp-shell Hamiltonians. Numerous detailed spectroscopic 
calculations have been reported. For example, in Ref. [|[, using a shell model approach, Martinez-Pinedo, Zuker, 
Poves and Caurier have performed full fp-shell calculations for the A=47 and A=49 isotopes of Ca, Sc, Ti, V, Cr 
£h ■ and Mn. They employed the KB3 interaction [8] with phenomenological adjustments and they performed complete 
diagonalizations to obtain very good agreement with the experimental level schemes, transition rates and static 
moments. Extensive discussions of fp-shell effective Hamiltonians and nuclear properties can be found in recent shell 
^ ■ model review articles [1, H, H, @] • 

Our interest in these nuclei stems from our goal to extend the ab-initio no-core shell model (NCSM) applications to 
heavier systems than previously investigated. Until recently, the NCSM, which treats all nucleons on an equal footing, 
had been limited to nuclei up through A = 16. However, in a recent paper d] we reported the first NCSM results 
for 48 Ca, 48 Sc and 48 Ti isotopes, with derived and phenomenological two-body Hamiltonians. These three nuclei are 
involved in double-beta decay of 48 Ca, and the interest in developing nuclear structure models for describing them 
is also related to the need for accurate calculations of the nuclear matrix elements involved in this decay. Our first 
goals were to see the limitations of such an approach applied to heavier systems and how much improvement one can 
obtain by adding phenomenological two-body terms involving all nucleons. In brief, the results were the following 0]: 
i) one finds that the charge dependence of the bulk binding energy of eight A=48 nuclei is reasonably described with 
an effective Hamiltonian derived from CD-Bonn interaction (Tlj. while there is an overall underbinding by about 0.4 
MeV/nucleon; ii) the resulting spectra are too compressed compared with experiment; iii) when isospin-dependent 
central terms plus a tensor interaction are added to the Hamiltonian, one achieves accurate total binding energies for 
eight A=48 nuclei and reasonable low-lying spectra for the three nuclei involved in double-beta decay. Only five input 
data were used to determine the phenomenological terms - the total binding of 48 Ca, 48 Sc, and 48 Ti along with the 
lowest positive and negative parity excitations of 48 Ca. 

In the present paper we extend our previous approach to the odd-A isotopes 47 Ca, 49 Ca, 47 Sc and 47 K, which differ 
by one nucleon from 48 Ca. One of our goals is to test whether the same modified effective Hamiltonian used for A=48 
isotopes, is able to describe these odd-A nuclei. A particular feature of the spectroscopy of these odd nuclei is that 
the spin-orbit splitting gives rise to a sizable energy gap in the fp-shell between the /y/2 and other orbitals (pi/2, P3/2, 
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/s/a) an d we wanted to see if this feature is reproduced in the NCSM where we have no input single-particle energies. 
Also, in spite of the differences in frameworks with and without a core, we wanted to compare our initial and modified 
Hamiltonian with a recent fp-shell interaction, the GXPF1, developed by Honma, Otsuka, Brown and Mizusaki [12 ]. 
We feel it is valuable to compare various fp-shell interactions in order to understand better their shortcomings and 
their regimes of applicability. From the comparison we present here, notable differences are evident. For example, 
our diagonal two-body T=0 matrix elements are more attractive and, while our initial and modified NN Hamiltonian 
fp-shell matrix elements are strongly correlated, there is much less correlation with the GXPF1 matrix elements, ft 
is worth mentioning that our interaction and Honma et al. GXPF1 interactions were also tested recently within the 
framework of spectral distribution theory in Ref. [l3j and sizable differences were demonstrated. 

Our paper is organized as follows: in Section 2 we give a short review of the NCSM approach and we refer the reader 
to the bibliography for more details. Section 3 is devoted to the presentation of our results. Binding energies, excitation 
spectra, single-particle characteristics, monopolc matrix elements and matrix element correlations are discussed in 
subsections along with corresponding figures. In the last Section we present the conclusion and the outlook of our 
work. 



II. NO-CORE SHELL MODEL 

The NCSM [3, Ei EE G3, EE Ell is based on an effective Hamiltonian derived from realistic "bare" interactions 
and acting within a finite Hilbert space. All A-nucleons are treated on an equal footing. The approach is both 
computationally tractable and demonstrably conver gent to the exact result of the full (infinite) Hilbert space. 

Initial investigations used two-body interactions [l4[ based on a G-matrix approach. Later, we implemented a 
similarity transformation procedure based on Okubo's pioneering work [20| to derive two-body and three-body effective 
interactions from realistic NN and NNN interactions. 

Diagonalization and the evaluation of observables from effective operators created with the same transformations 
are carried out on high-performance parallel computers. 

A. Effective Hamiltonian 

For pedagogical purposes, we outline the ab initio NCSM approach with NN interactions alone and point the reader 
to the literature for the extensions to include NNN interactions. We begin with the purely intrinsic Hamiltonian for 
the A-nucleon system, i.e., 

^ = T rel + V = if;fc|^+ £ Vvin-rj), (1) 

i<j i<j=l 

where m is the nucleon mass and Vn —rj), the NN interaction, with both strong and electromagnetic components. 
Note the absence of a phenomenological single-particle potential. We may use either coordinate-space NN potentials, 
such as the Argonne potentials 21] or momentum-space dependent NN potentials, such as the CD-Bonn [ill ]. 

Next, we add to (fT]) the center-of-mass Harmonic Oscillator (HO) Hamiltonian Hcm — Tqm + Ucm, where Uqm = 

convergence, the added -Hcm term has no influence on the intrinsic properties. 
However, when we introduce our cluster approximation below, the added Hcm term facilitates convergence to exact 
results with increasing basis size. The modified Hamiltonian, with pseudo-dependence on the HO frequency Q, can 
be cast as: 

A 

H A 2 = Ha + Hqm — ^ 

i=l 

Next, we introduce a unitary transformation, which is designed to accommodate the short-range two-body correla- 
tions in a nucleus, by choosing an anti-hermitian operator S, acting only on intrinsic coordinates, such that 

TL = e~ s H%e s . (3) 

In our approach, S is determined by the requirements that TL and H^ have the same symmetries and eigenspectra over 
the subspace K, of the full Hilbert space. In general, both S and the transformed Hamiltonian are A-body operators. 
Our simplest, non-trivial approximation to TL is to develop a two-body (a = 2) effective Hamiltonian, where the 
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upper bound of the summations "A" is replaced by "a", but the coefficients remain unchanged. We then have an 
approximation at a fixed level of clustering, a, with a < A. 



A ,A\ 



W = + 73^ £ VW., o) (4) 

i=l \a)\2! i 1 <i 2 <...<i a 



with 



V 12 ... a = e- sM H^e sM -$>, (5) 

i=l 

and S'^ a - ) is an a-body operator; = hi + h% + + . . . + h a + 14, and 14 = 53i<j ^i- We adopt the HO basis states 

that are eigenstates of the one-body Hamiltonian X^iLi hi- 

The full Hilbert space is divided into a finite model space ("P-space") and a complementary infinite space ("Q- 
space"), using the projectors P and Q with P + Q = 1. We determine the transformation operator S a from the 

decoupling condition Q a e~ s ° H^e s ° P a = and the simultaneous restrictions P a S^P a = Q a S^Q a = 0. The 
a-nucleon-state projectors (P a , Q a ) follow from the definitions of the A-nucleon projectors P, Q. 

In the limit a — > A, we obtain the exact solutions for dp states of the full problem for any finite basis space, with 
flexibility for choice of physical states subject to certain conditions (22J. This approach has a significant residual 
freedom through an arbitrary residual P Q -space unitary transformation that leaves the et-cluster properties invari- 
ant. Of course, the A-body results obtained with the a-body cluster approximation are not invariant under this 
residual transformation. An effort is underway to exploit this residual freedom to accelerate convergence in practical 
applications. 

The model space, P%, is defined by A m via the maximal number of allowed HO quanta of the A-nucleon basis 
states, Nm, where the sum of the nucleons' 2n + I < N m + A spsm i n = Am, and where A spsm i n denotes the minimal 
possible HO quanta of the spectators, nucleons not involved in the interaction. For example, 10 B, -/V spsm i n = 4 as 
there are 6 nucleons in the Op-shell in the lowest HO configuration and, e.g., A m = 2 + A max , where A max represents 
the maximum HO quanta of the many-body excitation above the unperturbed ground-state configuration. For 10 B, 
Am = 12, A m = 8 for an A max = 6 or "SMI" calculation. With our cluster approximation, a dependence of our 
results on A max (or equivalently, on A m or Am) and on Q arises. The residual A max and ft dependences will infer 
the uncertainty in our results arising from effects associated with increasing a and/or effects with increasing A max . 
In the present work, we retain the A max = basis space and hfl = lOMeV employed in Ref. @. 

At this stage we also add the term Hcm again with a large positive coefficient (constrained via Lagrange multiplier) 
to separate the physically interesting states with 0s CM motion from those with excited CM motion. We diagonalizc 
the effective Hamiltonian with the m-scheme Lanczos method to obtain the P-space eigenvalues and eigenvectors 
[ici| . All observables are then evaluated free of CM motion effects [l(| ■ In principle, all observables require the same 
transformation as implemented for the Hamiltonian. We obtain small renormalization effects on long range operators 
such as the rms radius operator and the B{E2) operator when we transform them to P-space effective operators at 
the a = 2 cluster level Il8l, 12-41 ] . On the other hand, when a=2, substantial renormalization was observed for the 
kinetic energy operator [23l|. and for higher momentum transfer observables [24| . 

Recent applications include: 

(a) spectra and transition rates in p-shell nuclei [25[ ; 

(b) comparisons between NCSM and Hartree-Fock [26[ ; 



(c) di- neutron correlations in the 6 He halo nucleus 271 ]; 

(d) neutrino cross sections on 12 C [281 ]; 

(e) novel NN interactions using inverse scattering theory plus NCSM (29l. l30l. [3l| ; 

(f) plus others in nuclear theory and quantum field theory (32j. 

We close this theory overview by referring to the added phcnomcnological NN interaction terms found adequate 
for obtaining good descriptions of A = 48 nuclei Q . Three terms are added - central modified gaussians with isospin 
dependent strengths and a tensor term. NCSM results obtained below with the modified Hamiltonian are referred 
to with "CD-Bonn + 3 terms" results. It is our hope that these terms accommodate, to a large extent, the missing 
many-body forces, both real and effective. This hope will be tested in the future when increasing computational 
resources allow larger basis spaces, improved a = 3 and a = 4 calculations as well as the introduction of true NNN 
and NNNN potentials. 
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III. RESULTS AND DISCUSSION 



A. Binding energies 



First we present the calculated total interaction energies (Hamiltonian ground state eigenvalues) in Fig. 1 which 
we compare with experiment. One observes that ground states calculated with our derived ab — initio H e ff lie above 
the experimental values by approximately 20MeV. This shift is similar to that observed in the case of all A = 48 
isotopes We note that with CD-Bonn we have nearly the same increase in binding from 47 Ca to 48 Cet as from 
48 Ca to 49 Ca, which signals a lack of subshcll closure. 

For the modified Hamiltonian (CD-Bonn + 3 terms) the NCSM produces reasonable agreement with experiment 
with deviations much less than 1% as seen in Fig. 1. There is a simple spreading of the theoretical ground states 
relative to experiment. In particular, we now observe the desired subshell closure condition where the increased 
binding from 47 Ca to 48 Ca significantly exceeds that from 48 Cet to 49 Ca. 
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FIG. 1: The experimental and theoretical ground state 
energy levels for A=47-49. The results in the second and 
third columns are labelled by their Hamiltonians. 



FIG. 2: Experimental and theoretical excitation energy lev- 
els for 49 Ca. Both CD-Bonn and CD-Bonn+3terms resultss 
are presented. 




FIG. 3: Experimental and theoretical excitation en- 
ergy levels for 47 Ca. Both CD-Bonn and CD- 
Bonn+3terms results are presented. 



FIG. 4: Experimental and theoretical excitation 
energy levels for 49 Sc. Both CD-Bonn and CD- 
Bonn+3terms results are presented. 
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FIG. 5: Experimental and theoretical excitation 
energy levels for 47 K. Both CD-Bonn and CD- 
Bonn+3terms results are presented. 
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FIG. 7: (color online) Comparison between CD- 
Bonn and GXPF1 of the monopole matrix elements 
V(ab;T)(A=48), which are shown by circles and squares, 
respectively. Lines are drawn to guide the eyes. The or- 
bit pair label "f7p3" stands for a = f 7 / 2 and b = P3/2, 
for example. 



FIG. 6: Experimental and theoretical levels that are 
dominantly single-proton and single-neutron particles 
or holes of 48 Ca. The levels are labeled by (n, 1, 2j), 
and the dashed lines are the Fermi energies. 
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FIG. 8: (color online) Comparison between CD- 
Bonn+3terms and GXPF1 of the monopole matrix 
elements V(ab;T)(A=48), which are shown by circles 
and squares, respectively. See the caption to Fig. 7. 



B. Excitation energy spectra 

The excitation energy spectra for 49 Ca, 47 Ca, A9 Sc and i7 K are shown in Figs. 2-5 respectively. In every case the 
ab — initio NCSM results with CD-Bonn arc far too compressed relative to experiment - a feature also seen in the 
A = 48 results [|. Here, we trace this primary defect to the inferred properties of the neutron orbits. That is, the 
incorrect ground state spin seen in Fig. 2 and the absence of a significant excitation energy gap in Fig. 3 indicate 
the spin-orbit splitting of the neutrons is insufficient to provide proper subshell closure at the neutron O/7/2 orbit. 
This defect is rectified in the results with CD-Bonn + 3 terms as seen by the corresponding spectra in Figs. 2 and 3. 
Similar tendencies have been seen before with valence G-matrix interactions and identified as a problem with the L 2 
dependence of the single-particle states 0, 0] • 

The CD-Bonn results in Figs. 4 and 5 are more difficult to interpret due to the glaring deficiencies just mentioned 
for the neutrons with the CD-Bonn Hamiltonian. We will show below that the proton shell closure is better established 
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FIG. 9: (color online) Comparison of the monopole 
matrix elements V(ab;T)(A=48) between CD-Bonn and 
GXPF1 (shifted to have an overall average monopole 
the same as CD-Bonn), which are shown by circles and 
squares, respectively. See the caption to Fig. 7. 
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FIG. 10: (color online) Comparison of the monopole ma- 
trix elements V(ab;T)(A=48) between CD-Bonn+3terms 
and GXPF1 (shifted to have an overall average monopole 
the same as CD-Bonn + 3 terms), which are shown by 
circles and squares, respectively. See the caption to Fig. 
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FIG. 11: (color online) Correlation of V(abcd; JT) 
(A=48) matrix elements between CD-Bonn+3terms and 
CD-Bonn. The matrix elements of T=0 and T=l are 
shown by open circles and crosses, respectively. There 
are no monopole shifts. The straight line represents a 
linear fit to all the matrix elements. 



FIG. 12: (color online) Correlation of V(abcd; 
JT)(A=48) between CD-Bonn and G. See the cap- 
tion to Fig. 11. 



with CD-Bonn. This supports the assertion that the main deficiencies seen in the third columns of Figs. 4 and 5 are 
indeed likely to reside with the inferred neutron spin-orbit splitting problem. 

The modified Hamiltonian provides greatly improved spectra for all four nuclei as seen in the second columns 
of Figs. 2-5. It is to be noted that these nuclei were not involved in the fitting procedure used to determine the 
parameters of the added phenomenological terms. Perhaps the most significant remaining deficiency is the incorrect 
ground state spin for 47 K as seen in Fig. 5. This is the first case of a nucleus in the region of A — 47 to A — 49 (12 
nuclei studied to date) where we did not obtain the correct ground state spin with CD-Bonn + 3 terms Hamiltonian. 



C. Single-particle characteristics 



In order to better understand the underlying physics of our NCSM results, we investigated the single-particle-like 
properties of our solutions. 
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FIG. 13: (color online) Correlation ol the ma- 
trix elements V(abcd; JT)(A=48) between CD- 
Bonn+3terms and GXPF1. See the caption to Fig. 
11. 




FIG. 14: (color online) Correlation of the ma- 
trix elements V(abcd; JT)(A=48) between CD- 
Bonn+3terms and GXPF1, where we retain only the 
off-diagonal matrix elements, i.e. the 56 matrix el- 
ements that cannot contribute to a single-particle 
Hamiltonian (see text). See the caption to Fig. 11. 



In a simple closed-shell nucleus, we expect the leading configuration of the ground state solution in our m-scheme 
treatment to be a single Slater determinant. Single-particle (or hole) excitations should be easily identified by the 
character of their leading configurations - i.e. a single-particle creation (or destruction) operator acting on the ground 
state Slater determinant of the reference nucleus 48 Ca. For our odd-mass nuclei, this is the character we seek. That 
is, we take the standard phenomenological shell model configuration of a single Slater determinant with a closed 
sd-shell for the protons and a closed f-j/2 subshell for the neutrons and look for the appropriate states which have 
a single nucleon added to (or subtracted from) that Slater determinant. We accept states as "single-particle-like" 
when we find one with a leading configuration having more than 50% probability to be in the simple configuration 
just described. When the majority weight is distributed over a few states, we use the centroid and we discuss those 
cases in some detail below. We were not successful in locating all the expected single-particle-like and single-hole-like 
states. That is, those absent from our presentation below were spread among a large number of eigenstates. 

For a closed-shell nucleus (Z,N) the single-particle energies (SPE) for states above the Fermi surface are related to 
the binding energy differences: 



= BE(Z, N) - BE*(Z + 1,N), 
= BE(Z, N) - BE*(Z, N + 1). 
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The SPE for sates below the Fermi surface are given by 
e< = BE*(Z - 1,N) - BE(Z, N), 
and 

e< = BE*{Z, N-l)- BE{Z, N). 

The BE are ground state binding energies which are taken as positive values, and e will be negative for bound states. 
(BE* = BE — E x ) is the ground state binding energy minus the excitation energy of the excited states associated 
with the single-particle states. 

Experimental SPE's and the results of our analysis are shown in Fig. 6. The experimental SPE's for protons and 
neutrons follow B.A.Brown's analysis [4(. To guide the eye, we draw a horizontal line to indicate the vicinity of the 
Fermi surfaces for the protons and neutrons. 

Fig. 6 shows that proton shell closure is established with both Hamiltonians, the CD-Bonn and the CD-Bonn+3 
terms. The correct energy locations arc better approximated with the modified Hamiltonian. Fig. 6 also shows that 
neutron subshell closure only appears with modified Hamiltonian. Here, the ordering is correct but the states are 
considerably more spread out compared with experiment. 

Let us consider some of the details underlying the single-particle-like states. The situation for the 1^3/2 or "lp3" 
state in the left panel of Fig. 6, the proton single-particle state in 49 Sc, with the modified Hamiltonian, is quite 
interesting. It appears that this state is mixed over several excited states in the spectrum. We can take the strength 
spread over several states and construct a centroid for this lp3 state by a weighted average over the states carrying 
that strength. Here are the relevant input ingredients. 

The first excited state of 49 Sc is a 3/2~, as seen in the second column of Fig. 4, with about 51% of the occupancy 
of the lp3 state. Its eigenvalue is -425.151 MeV compared to a ground state of -428.365 MeV. The 18th state in the 
49 5c spectrum is also a 3/2~ with 28% of the occupancy of the lp3 state. Its eigenvalue is -422.803MeV. The 24th 
state is also a 3/2~ with 21% of the occupancy of the lp3 state. Its eigenvalue is -422.440MeV. 

Thus, to a good approximation, the lp3 strength is spread over these three states. We will identify the weighted 
average [0.51 x (-428.365) + 0.28 x (-422.803) + 0.21 x (-422.440)] = -423.79 as the centroid of the single particle 
lp3 state which we then include accordingly in the second column of the figure. 

For the proton hole states with the modified Hamiltonian, we perform a detailed search up to excitation energies of 
about 14 MeV in the 47 K spectra. It appears that the 0^5/2 single-hole state is spread among many states with the 
largest observed concentration on the 5/2 + state at -386.17McV (13.36 MeV of excitation energy). Here, we find a 
single J 71 " = 5/2 + state in 47 K with 30 % 0g?5/2 vacancy and we assign this state to our 0^5/2 single-hole state. Most 
of the 0c?5/2 strength, however, was not observed among the limited number of converged eigenstates. 

Let us consider the 49 Ca results with the modified Hamiltonian in the upper right panel of Fig. 6. The ground state 
is approximately a pure [(lp3/2) 1 (0/7/2) 8 ] configuration. We note that the spacing for the subshell closure is in good 
agreement with experiment while there is a shift of a couple MeV towards more binding in the model as previously 
indicated in Fig. 1. A nearly pure lpi/2 single-particle state is obtained at 5.235 MeV excitation energy and an extra 
low-lying 7/2~ appears with 2p-lh character (see Fig. 2). Our lowest-lying 5/2~ consists of 2p-lh character relative 
to subshell closure. 

We contrast the modified Hamiltonian's results for the 49 C 'a ground state with those obtained using the ab — initio 
CD-Bonn where [(lp3/2) 4 (lpi/2) 2 (0,/V/2) 3 ]i/2- is the dominant configuration reflecting again the inadequacies of the 
neutron single-particle properties. 



D. Monopole matrix elements V(ab;T) 

The monopole matrix element is defined by an angular momentum average of coupled doubly-reduced two-body 
matrix elements: 

For our NCSM Hamiltonians the " V" appearing in Eqn. [6] signifies the full 2-body intrinsic-coordinate Hamiltonian, 
T r ei + Veff, except that we omit the Coulomb interaction from this analysis. 

We examine the monopole character of our initial CD-Bonn Hamiltonian and we note some similarities and differ- 
ences from the GXPF1 interaction[l2j as shown in Fig. 7. The see-saw shapes of the two Hamiltonians in Fig. 7 
are similar but our Hamiltonian is shifted towards less attraction. Given the many differences between the respective 
theoretical starting points, the ab — initioH e ff for the NCSM and the G-matrix for GXPF1 the different bare NN 
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interactions, etc., the similarities observed in Fig. 7 are remarkable. In order to summarize a comparison of the 
underlying theoretical interactions, we list in Table la simplified overview of their differences and similarities. 

For a more detailed comparison of the interactions, we present the fp-shell matrix elements applicable to the 
present investigation in Tables 2 and 3. For convenience in finding the major differences, we present two columns 
of key differences in the matrix elements: "diffl" represents the difference between our ab — initio H e ff and our 
modified H e ff, and, "diff2" represents the difference between our modified H e ff and the GXPF1 interaction. While 
diffl shows magnitudes that only occasionally exceed 1 MeV, diff2 shows magnitudes approaching 2.6 MeV. This type 
of comparison suggests that our solution for the modified Hamiltonian remains closer to our initial H e f f derived from 
CD-Bonn, than it is to the fitted Hamiltonian, GXPF1. 

Fig. 8 presents a similar comparison of the monopole character in the fp-shell of the two phenomenological Hamil- 
tonians, GXPF1 and CD-Bonn + 3 terms. Overall, the changes in the monopole character due to the addition of the 
phenomenological terms to our i? e // of Fig. 7 appear somewhat larger for the T=l monopole than for T=0. The 
effect of "+3 terms" is to increase the T=0 and T=l splitting of six of the monopoles while the 4 remaining T=0 and 
T=l monopole splittings are reduced. 

In order to better visualize the similarities of the fp-shell matrix elements, we present in Figs. 9 and 10 the same 
comparisons shown in Figs. 7 and 8, respectively, with an overall shift of the GXPF1 monopole matrix elements so 
that the average over all monopole matrix elements is the same for the two Hamiltonians. Specifically the average 
shift of T=0 and T=l monopoles for GXPF1 in Fig. 9 is 0.899768 MeV, while for GXPF1 in Fig. 10 it is 0.83485 
MeV. It is now evident that both the CD-Bonn and CD-Bonn + 3 terms have monopoles with less T — and T = 1 
splittings. 

E. Matrix element correlations 

We present in Figs. 11-14 the correlations between pairs of fp-shell interaction matrix element sets. With Fig. 11, 
we observe the high degree of correlation between the 195 matrix elements of our starting Hamiltonian, CD-Bonn, and 
our modified Hamiltonian, CD-Bonn + 3 terms. This indicates that, for the most part, our Hamiltonian is minimally 
modified by the addition of the phenomenological terms. Such a high correlation is reminiscent of the high correlations 
seen between GXPF1 and its starting interaction, the G-matrix [12j ]. 

It is then very interesting to observe in Fig. 12 the lack of correlation between our starting Hamiltonian, CD-Bonn, 
and the G-matrix underlying the GXPF1 interaction. This lack of correlation reflects the major differences in the 
underlying theories that are summarized in Table 1. 

Furthermore, we find minimal correlation between CD-Bonn + 3 terms and the full GXPF1 as seen in Fig. 13. 
This indicates the likely sensitivity to the starting Hamiltonians in the fitting procedures and to the differences in the 
NCSM compared to a valence shell model approach. 

Finally, in order to isolate the true off-diagonal 2-body interaction effects from those that may lead to contributions 
to the single-particle energies, we eliminate some of the matrix elements from the comparison. In Fig. 14 we 
present the correlation of matrix elements V(abcd; JT)(A=48) between CD-Bonn+3terms and GXPF1, where we 
retain only those that cannot contribute to a single-particle Hamiltonian in leading order. That is, we eliminate all 
two-body matrix elements where at least one single-particle-state (sps) of the bra equals a sps of the ket. There 
are 56 remaining two-body matrix elements. Clearly, the correlation improves. This suggests that those eliminated 
matrix elements would generate single-particle properties, in leading order, different from the single-particle properties 
embodied in GXPF1 plus its associated single-particle energies. Comparisons of spectra and other properties with 
these Hamiltonians as one proceeds further from A=48 could shed more light on their differences. 

IV. CONCLUSIONS AND OUTLOOK 

We have presented an initial NCSM investigation of the spectral properties of the A = 47 and A = 49 nuclei 
that are one nucleon away from doubly-magic 48 Ca. We have shown that the NCSM with a previously introduced 
modified Hamiltonian produces spectral properties in reasonable accord with experiment. Shell closure properties 
are obtained and a path has been opened for multi-shell investigations of these nuclei within the NCSM. We are 
undertaking such additional investigations. Also, for a better understanding of various fp-shell interactions we made 
a comparison between our initial and modified fp-shell matrix elements in the harmonic oscillator basis with the 
GXPF1 interaction [12| and we found notable differences. Our initial and modified NCSM H e ff matrix elements in 
the fp-shell are strongly correlated. However, the same matrix elements of the modified NCSM Hamiltonian appear to 
lack a significant correlation with the GXPF1 matrix elements (see Fig. 13). We found evidence (Fig. 14) suggesting 
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that significant differences in single-particle properties may underly some of the distinctions between our H e f f and 
the GXPF1 interaction. Additional applications could reveal those distinctions in greater detail. 
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Reguires perturbative corrections 






and raises known convergence issues 
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Starting energy dependence 
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Single-particle spectra dependence 
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No 


A-dependence 
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TABLE I: Overview of the differences and similarities of the two theoretical approaches that underlie the Hamiltonians whose 
matrix elements are compared in this work. 
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TABLE II: Comparison of selected two-body matrix elements V(abcd; JT) (Mev) (A=48) for which the difference between our 
interaction is large. Diffl represents the difference between GXPF1 and G and diff2 is the difference between CD-Bonn+3terms 
and CD-Bonn. 
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TABLE III: Comparison of the fp-shell two-body matrix elements V(abcd; JT) (Mev) (A=48) 
interaction GXPF1 is taken from Ref. Diffl represents the difference between CD-Bonn+3t' 

is the difference between CD-Bonn+3terms and GXPF1. 
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